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i — i Abstract 

We study thermal leptogenesis from decays of the electroweak 
triplets that mediate neutrino masses in type-II and type-Ill see- 
saw. We find that Sommerfeld corrections reduce the baryon 
asymmetry by ~ 30%, and that successful leptogenesis needs 
cn triplets heavier than 1.6 TeV, beyond the discovery reach of LHC. 

1 Introduction 

o 

Leptogenesis might be produced by decays of scalar (type-II see-saw) or fermion (type-Ill 
see-saw [1]) SU(2)x triplets with mass M. The Boltzmann equation for their abundance 
contains their annihilation rate 7^4 (number of annihilations per spacetime volume at tem- 
perature T): the final baryon asymmetry depends on the value of 7,4 in the non-relativistic 
limit, at T ~ M/ln(Mp\/M). Non-relativistic scatterings among gauge-charged particles 
are affected by non-perturbative corrections: the states involved in the annihilation no 
longer are plane waves when the kinetic energy is comparable to the electroweak potential 
energy [2]. In this paper we study how such Sommerfeld corrections affect the final baryon 
asymmetry. 

In section 2 we use group theory to derive a simple formula for non-abelian Sommerfeld 
corrections. In section 3 we apply it to type-Ill see-saw and in section 4 to type-II see-saw. 



1 



In both cases we also extend computations of thermal leptogenesis down to values of 
M ~ TeV accessible to the LHC collider, motivated by recent papers that studied how 
type-II or III see-saw can be searched for at LHC. However, we find a model-independent 
lower bound on M that limits this possibility. In section 5 we discuss our results and 
conclude. 



2 Sommerfeld corrections 

We denote as T the scalar or fermion Triplet that decays producing the lepton asymmetry. 
Following [3, 4] the final baryon asymmetry is parameterized as ris/n^ = — 0.029i]e, where 
the efficiency 77 encodes the dynamics and e the amount of CP violation. The Boltzmann 
equation that dictates the evolution of the triplet abundance contains the T decay rate 
7d and the TT* annihilation rate 7^: 

TT dY ( Y \ ( Y 2 \ 

where z = M/T, H is the Hubble rate at temperature T; Y = nq-/s for the real fermion 
triplet and Y = (riq- + nq-*)/s for the complex scalar triplet; Y eq is the value that Y would 
have in thermal equilibrium; n is the number density; s = 2n 2 g* s T 2 ' /45 with g* s = 106.75 
is the entropy density of SM particles; 7^ is the thermal average of the annihilation cross 
section summed over all initial- and final-state Lorentz and gauge indices: 



T 



r ds ^W^Ws) where a A = [ dt^ 1 -^, (2) 
Jam 2 \ 1 J J ^ Sits 



and Ki is a Bessel function; 

Explicit expressions for a a are given in appendix A of [5] for a generic SU(2)i<g)U(l)y 
multiplet. Similarly to the case of Cold Dark Matter freeze-out, the value of 7^ is relevant 
in the non relativistic limit, T/M w 1/lnMpi/M 1, so that a a can be approximated 
as the sum of s-wave, p-wave, etc contributions: 

a A = c s (3 + c p (3 3 + 0((3 5 ) (3) 

where (3 = \Jl — AM 2 / s is the triplet velocity in the TT* center of mass frame. The 
corresponding annihilation rate is 



MT 3 e -2M/T 
327T 

For the fermion triplet one has 



~r<J 



(4) 



lllflrt 51^ _ 

V' ^=^- (5) 
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For the scalar triplet 7^ is well approximate by just its s-wave coefficient 

9gj + 12gfg| + 34 

(hypercharge is normalized such that Y = 1/2 over lepton doublets). 



(6) 



We now compute the non-perturbative electroweak Sommerfeld corrections to c s . Scat- 
terings among charged particles are distorted by the Coulomb force, when their kinetic 
energy is low enough that the electrostatic potential energy is relevant. This leads e.g. 
to significant enhancements of the annihilation cross section (attractive force) or 

to significant suppressions of various nuclear processes (repulsive force). The Sommerfeld 
correction to s-wave point-like annihilations (e.g. \T n + — > 77) can be computed as [2] 
S = \ip(oo)/^(0)\ 2 , where ip{r) is the (reduced) s-wave-function for the two-body state 
with kinetic energy K = M/3 2 , that in the non-relativistic limit satisfies the Schr6dine;er 
equation 

with outgoing boundary condition ^//(oo)/^>(oo) ~ iM/3. For a single abelian massless 
vector with potential V = a/r the result is a = SVperturbative where the Sommerfeld 
correction is [2] 

S(x) = -f^- x= a o . (8) 

Here a < describes an attractive potential that leads to an enhancement S > 1, and 
a > describes a repulsive potential that leads to S < 1. We will need the thermally 
averaged value of the Sommerfeld correction: 

J P 2 e- M ^/ T S(a/P)dp = a_ 

where (3t = \Jt/M is the characteristic velocity in the thermal bath and the upper 
integration limit on (3 is any value much larger that (3t- The function St is computed 
numerically, with the qualitative result St(x) ~ S(x). 

The generalization to non-abelian massive vectors was discussed in [6, 5] and needs a 
long list of potential and annihilation matrices, such that the matrix Schrodinger equa- 
tion (7) can be solved only numerically. We here show how in the SU(2)£,-invariant limit 
M ^> Mw t z the Sommerfeld correction can be analytically computed as a sum of abelian- 
like cases. For any simple gauge group with coupling a = g 2 /4:ii, the potential between 
two particles in representations R and R' at distance r is 

V=-Y.T a R®T a R , (10) 
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where a runs over the adjoint, T R are the generators in the representation R, and the 
tensor product (g) indicates that V has 4 gauge indices. Unlike in [5] we do not include 
(anti)symmetrizations and their corresponding normalization of the 2-body states: their 
only effect is enforcing V — for the 2-fermion states forbidden by the Pauli exclusion 
principle, which can be less formally imposed by hand. Group theory allows to decompose 
the tensor product of the representations R and R' into a sum of irreducible representa- 
tions Q: 

R®R> = @Q, T° m , ee T% ® l R , + l R ®T a RI =Y. T% (11) 
Q Q 

so that, recalling the definition T R -T R = C R 1 R of the quadratic Casimir C R , gives: 
V = ~ Y.^R'-T a mK -T a -T a ®l RI -l R ®T§.T^) = ^ (£C Q l Q -C R l-C R d) (12) 

A Q 

where 1q is the projector along the subspace Q, so that the matrix V is diagonal along 
each irreducible representation Q. The Casimir of the SU(2) irreducible representations 
with dimension n is C n = (n 2 — l)/4, so that a two-body state n © n with total isospin 
iV < 2n — 1 has potential V = (N 2 + 1 — 2n 2 )ct2/8r. The two-body state of the (scalar 
or fermion) triplets involved in leptogenesis decompose as 3 ® 3 = lg © 3^ © 5s, and the 
potentials are V = —2a 2 /r for the singlet state, V = —a^jr for the triplet state, and 
V = a 2 /r for the quintuplet state. 

Finally, one might worry that the Sommerfeld correction, computed for scattering in 
vacuum, does not apply for annihilations in the cosmological plasma. This is not the 
case because is only needed when the TT* annihilation rate becomes smaller than 
the expansion rate H ~ T 2 /Mp\. in view of the Planck suppression H and thereby 7a 
is so small that annihilations are rare enough that the vacuum approximation holds. 
Furthermore, the T life-time r ~ 1/(A 2 M) is longer enough than the Coulomb time-scale 
~ 1/(M(3 2 ) provided that A -C /3, where A denotes the small coupling(s) present in type- 
Ill or type-II see-saw. Finally, thermal effects generate a Debye mass for the vectors 
(m 2 = 1 1(72^ 2 / 6 for the SU(2)l vectors), screening the Coulomb potential into a Yukawa 
potential, e~ mr /r: the results of [5] show that the Debye mass m is small enough that we 
can neglect it, approximating vectors as massless. 

3 Type-Ill see-saw: leptogenesis from a fermion triplet 

The two body TT states can be classified according to their quantum numbers (/, 5*, L), 
where / = {1,3,5} is the total isospin, S = {0, 1} is the total spin, L is total orbital 
angular momentum. Restricting to the dominant s-wave annihilations, L = 0, the states 
allowed by quantum statistics are (I,S) = (1,0), (3,1) and (5,0). Their annihilations 
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Figure 1: Annihilation rate 7^ for the fermion and scalar triplet for M = 10 5 GeV. The 
'tree' line shows the tree-level result, the dashed line shows the s-wave approximation, 
and the 'Sommerfeld' line includes Sommerfeld corrections. Leptogenesis is dominantly 
produced at z = M/T ~ 20. 



into couples of SM vectors (W a and Y of SU(2)^ and hypercharge), of SM fermions ^ 
and of SM Higgs doublets H are restricted by their (J, S, L) quantum numbers. Taking 
into account the potentials computed in section 2, the total Sommerfeld-corrected s-wave 
annihilation rate is: 1 

c s = ^is T (-2a 2 ^) + 5 As T ( a2 ^) + (1 + _L)^ T (-a 2 ^) (13) 

7T Z7T Z4 7T 
s v ' " v ' s v ' 

(I,S)=(lfl)^W a W a (I,S) = (5,0)^W a W b (J,S)=(3,l)-»Wf,M 

-'^We compare with previous results. Neglecting Sommerfeld corrections, St — > 1: i) eq. (13) agrees 
with eq. (29c) of [3] (where a factor 2 was however missed in the Boltzmann equation in front of ja the 
factor 1/24 coming from annihilations into HH* was neglected; the g\ factor was unproperly typed as g\ 
in the numerical code); ii) eq. (13) differs by order one factors from eq. (4.9) of [7]; iii) eq. (13) agrees with 
analogous computations of the relic abundance of wino-like dark matter. Some dark matter studies also 
included Sommerfeld corrections [6, 5], but a direct comparison is not immediate (dark matter studies 
included SU(2)i-breaking effects and thereby used electric charge instead of isospin to classify states, 
obtaining complex expressions); the numerical results for 7^ can be compared in the limit M 3> Mzjoii 
and agree. Some recent experimental results suggested the possibility of Sommerfeld enhancememnts 
from new gauge interactions in the Dark Matter sector [8]: our result in eq. (11) holds for any gauge 
group and allows to study such effects. 
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Figure 2: Iso-contours of the efficiency 77 of thermal leptogenesis from decays of a fermion 
triplet (left, type-Ill see-saw) and a scalar triplet (right, type-II see-saw, in the case of 
equal branching ratio into leptons and into higgses). The regions shaded in green are 
allowed by the model- dependent bound on the CP asymmetry generated by the neutrino 
mass operator (LH) 2 . 



where the text under each contribution to the rate indicates the corresponding annihilation 
processes. Eq. (13) differs from the perturbative result of eq. (5) by an order one factor 
at T<Mg 2 2 . 

Fig. la shows our result for the adimensional combination 7^e 2z z 3 /M 4 , that becomes 
constant in the non-relativistic limit z = M/T ^> 1, where s-wave annihilations domi- 
nate. The final baryon asymmetry roughly depends on 7^4 at z ~ 20: the Sommerfeld 
enhancement to the s-wave contribution is more important than including />wave and all 
other L 7^ annihilations. 

Fig. 2a shows the contour plot of the efficiency 77 as function of M and of mi, which, 
as usual, is the lightest triplet contribution to neutrino masses. Considering e.g. the 
point M = 10 5 GeV and rh\ = 10~ 5 eV, Sommerfeld corrections reduce the efficiency from 
2.6 1(T 8 to 2.0 10~ 8 . 

For large enough rh\ the decay rate 7^ is larger than the annihilation rate 'Ja in the 
relevant temperature range T ~ M/20, so that the efficiency 77 no longer depends on 7^4, 
and gets the value typical of type-I see-saw, largely independent on M. In this region 
Sommerfeld corrections, that enhance 7a, do not affect the final asymmetry. 
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Figure 3: Values of the CP asymmetry \e\\ < 1 needed to get the observed baryon asym- 
metry from decays of a fermion triplet (left) and a scalar triplet (right), as function of 
its mass M for the values of fhi superimposed to the curves. We see that successful 
leptogenesis needs M > 1.6 TeV. 



As in [3] we shaded in green the region where thermal leptogenesis can produce the 
observed baryon asymmetry, assuming that the CP asymmetry is generated by two other 
triplets so heavy that their effects are fully encoded in the dimension-5 operator (LH) 2 , 
constrained by neutrino masses. We see that the resulting bound on the CP asymmetry [3] 
implies the model-dependent lower bound M > 3 10 10 GeV. 

Various recent papers considered type-Ill and type-II see-saw with M so small that 
triplets can be probed by colliders: LHC experiments are going to probe the region 
100 GeV < M < 1 TeV. We here find that this range of M is not compatible with thermal 
leptogenesis, that needs M > 1.6 TeV because at lower M the efficiency becomes so small 
that the baryon asymmetry ub/u^ = — 0.0297/ei [3] remains smaller than the observed 
value, ne/n 7 ~ 6.15 10~ 10 [9], even if the CP asymmetry is maximal, |ei| = 1 (resonant 
leptogenesis [10] allows to realize an order one CP-asymmetry, although unity is not 
reached). This issue is better illustrated by fig. 3a, that also allows to see how the lower 
bound on M depends on rh\. 

Let us now discuss why the efficiency strongly decreases when M is below 3 TeV. The 
Higgs starts to acquire its vev via a first order phase transition at the critical temperature 
T cr pa 1.2rrih- Sommerfeld corrections have to be computed as in [6, 5], but this is a minor 
detail. The key issue is that SU(2)^-breaking suppresses the sphaleron rate 7 sp haieron [11], 
so that the lepton asymmetry ceases to be converted into the baryon asymmetry below 
the temperature T dec pa 80 GeV + 0A5m h [11]. The approximations for T cr and T dec hold 
for the light Higgs mass suggested by present data, 115 GeV < rrih ^ 200 GeV [12]; we here 
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assume = 120 GeV. In order to quantitatively include the sphaleron decoupling effect 
it is convenient to write the Boltzmann equations for the B — £ and the B asymmetries: 



(14a) 

c^0.52 (14b) 

The first equation is the standard one [3], as sphalerons conserve B — £? The second 
equation tells the final baryon asymmetry, and is approximatively solved by 

Y B {T < T dcc ) w (1 + l/c)F B _ £ (T dcc ). (15) 

The evolution of the various asymmetries is illustrated in fig. 4a, where we see that the 
observed baryon asymmetry is reproduced with a maximal CP asymmetry = 1 for 
M = 1.6 TeV and rhi = 0.06 eV, and it is dominantly produced at low temperature 
T ~ M/10, where neglected higher order processes are expected correct it by < 10%. 

Fig. 3a shows how the CP asymmetry E\ needed to get the observed baryon asymmetry 
depends on M for several values of rh\. Since the lepton asymmetry is dominantly gener- 
ated by triplet decays at T ~ M/20, the sphaleron decoupling effect becomes significant 
at M>20T dec ~ few TeV, preventing successful baryogenesis if M > 1.6 TeV. 



„ dY B -c ( y ,\ Yb-cId 



dY 

SHZ—— = -7sphaleron(Vf5C - (1 + c)Y B - C ) 



4 Type-II see-saw: leptogenesis from a scalar triplet 

The two body TT* states have total isospin I = {1, 3, 5}, total spin S = 0; and we focus 
on s-wave annihilations so that L = 0: at tree level such states can only annihilate into 
two SM vectors, that can have isospin / = {1, 3, 5}. 



c s = 2 ^ 4 + 3^/2 ^ ( _ (2a2 + + ^As T ((a 2 - a Y )^z) + 6 -^S T (-(a 2 + a Y )^z) 

7T ZTT 7T 
v v ' v v ' v v ' 

I=l^W a W a .YY I=5^W a W b I=3^W a Y 

(16) 

In the limit St — *■ 1 this expression for 7^ reduces to the one in [4]. 

Our results for the type-II see-saw are presented in fig. lb, 2b, 3b, that closely resemble 
the corresponding figures la, 2a, 3a already presented in the type-Ill section. 



including flavor details the first equation gets replaced by three equations for the partial B/3 — C e ,n,T 
asymmetries [13]. If fermion triplet(s) will be discovered at LHC, it should be possible to reconstruct the 
flavor structure of their couplings [14] and a more precise analysis will become worthwhile. 
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Figure 4: Evolution ofY L , Y B , Y — Y eq for M = 1.6 TeV, m x = m atm = 0.06 eV, \e x \ = 1. 



Fig. lb shows how much Sommerfeld corrections enhance the annihilation rate 7,4, 
thereby decreasing the efficiency. Fig. 2b shows the contour plot of the efficiency rj as 
function of M and of fh\, assuming equal T branching ratios into LL and into HH, 
Bl = Bh = 1/2. Considering e.g. the point M = 10 5 GeV and rh\ = 10~ 5 eV, Sommerfeld 
corrections reduce the efficiency from 2.9 10~ 8 to 2.1 10 -8 . 

As in [4] we shaded in green the region of fig. 2b where thermal leptogenesis can 
produce the observed baryon asymmetry, assuming that the CP asymmetry is generated 
by the dimension-5 operator (LH) 2 , constrained by observations of neutrino masses. 

Finally, we study the model-independent leptogenesis lower bound on M. The ef- 
ficiency is higher when Bl S> Bh or Bh 3> Bl [4]. However the model-independent 
unitary bound on the CP asymmetry in T decays, \el\ < 2 min(I?£, Bh) [4] has an oppo- 
site and stronger dependence on Bl, Bh, such that the lower bound on M is saturated 
for Bl = Bh = 1/2 (we are not aware of explicit models that can saturate this bound). 
As in the previous section, the key point is taking into account sphaleron decoupling, by 
adding the sphaleron eq. (14b) to the standard equations of [4]. The resulting bound 
M > 1.6 TeV is illustrated in fig. 3b, which is qualitatively and quantitatively similar to 
fig. 3a. 

Fig. 4b shows the evolution of the various asymmetries for M = 1.6 TeV, rh\ = 0.06 eV, 
and a maximal CP asymmetry \el\ = 1- 



9 



5 Conclusions 



Type-II and type-Ill see-saw generate neutrino masses from particles with electroweak 
gauge interactions: scalar or fermion SU(2)^ triplets. Previous works showed that these 
triplets can lead to successful thermal leptogenesis [3, 4, 7] as %/n 7 = — 0.029e^ = 
6 • 1CT 10 . The triplet annihilation rates (that keep triplet abundances close to thermal 
equilibrium apparently violating the out-of-equilibrium Sakharov condition for baryogen- 
esis) can be small enough to allow a large enough efficiency 77. 

In the present work we showed that Sommerfeld corrections, that account for long- 
range non-abelian electrostatic interactions, enhance the annihilation rate, decreasing the 
efficiency by ~ 30%. Fig. 2a and b show the efficiency as function of the lightest triplet 
mass M and of the lightest triplet contribution rhi to neutrino masses. 

The efficiency steeply decreases at M<3TeV because, when at T ~ M/20 triplet 
decays produce the lepton asymmetry, sphalerons no longer convert it into the baryon 
asymmetry. Indeed at T < the electroweak symmetry starts to be broken suppressing 
the sphaleron interaction rate. We find that baryogenesis via type-II or III thermal 
leptogenesis is only possible at M> 1.6 TeV. This bound on M becomes stronger if rh\ 
is smaller than 1 eV, or if the Higgs is heavier than its minimal value, rrih = 115 GeV. 

Previous studies [15, 14] found that see-saw triplets lighter than about 1 TeV give 
detectable effects at the LHC collider. It is therefore interesting to study if /how the 
absolute leptogenesis lower bound on M can be evaded. One possibility is a new source of 
baryogenesis at the weak scale, for example new physics that makes the electroweak phase 
transition of second order and provides a new source of CP violation [16], a scenario which 
is already strongly constrained but not yet excluded. The alternative possibility is another 
heavier baryogenesis mechanism: here the problem is that the baryon asymmetry gets 
washed out by triplets unless their L-violating interactions are slower than the expansion 
rate. In the fermion triplet case, this happens if fh\ <C 10~ 2 eV (which implies significantly 
displaced decay vertices at LHC [14]) or if the triplet is negligibly coupled to some lepton 
flavor. However, neutrino oscillations need neutrino masses heavier than 10~ 2 eV and with 
large mixings among flavors. Scalar triplets offer one more possibility [4]: L- violating 
rates are suppressed if triplet decay rates into LL or HH are much different, but this 
suppression would hold in cosmology as well as at LHC. 

We obtained the bound M > 1.6 TeV in type-II and type-Ill see-saw, and we expect 
a similar lower bound on M in any model where the particle responsible for leptogenesis 
has interactions large enough to be detectably produced at LHC. In particular, in type- 
I see-saw models where baryogenesis is produced by decays of a right handed neutrino 
charged under new massive vectors, one can maybe find a situation where the new vectors 
are light enough to lead to detectable signals at LHC, and heavy enough to suppress their 
contribution to [17]. 
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We conclude discussing the bound on M from the string-anthropic-multiverse point 
of view that recently received significant attention [18]. Assuming that baryogenesis is 
impossible in the vast majority of string models because they predict M <C (h), one 
could argue that it is more likely to find a M close to the anthropic bound that makes 
baryogenesis possible. If this possibility will be confirmed by future data, one might even 
argue that the weak scale is much below the Planck scale because it cannot be larger 
than M, which is generated exponentially below the Planck scale by string instanton 
effects [19]. 

Acknowledgements I thank Antonello Polosa and Marco Bochicchio for helpful dis- 
cussions about group theory. 

References 

[1] R. Foot, H. Lew, X.-G. He and G.C. Joshi, Z. Phys. C44 (1989) 441;. 

[2] A. Sommerfeld, Ann. Phys. 11 257 (1931). 

[3] T. Hambye et al., Nucl. Phys. B695 (2004) 169 [arXiv:hep-ph/0312203]. 

[4] T. Hambye, M. Raidal, A. Strumia, Phys. Lett. B632 (2006) 667 [arXiv:hep-ph/0510008]. 

[5] M. Cirelli, A. Strumia and M. Tamburini, Nucl. Phys. B787 (2007) 152 [arXiv:0706.4071]. 

[6] J. Hisano, S. Matsumoto, M. Nagai, O. Saito, M. Senami, Phys. Lett. B646 (2007) 34 
[arXiv:hep-ph/0610249]. 

[7] W. Fishier, R. Flauger, arXiv:0805.3000. 

[8] M. Cirelli, M. Kadastik, M. Raidal, A. Strumia, arXiv:0809.2409; N. Arkani-Hamed, D.P. 
Finkbeiner, T. Slatyer, N. Weiner, arXiv:0810.0713. 

[9] For recent reviews see: A. Strumia and F. Vissani, arXiv:hep-ph/0606054. M.C. Gonzalez- 
Garcia, M. Maltoni, arXiv:0704.1800. 

[10] M. Flanz, E.A. Paschos, U. Sarkar and J. Weiss, Phys. Lett. B389 (1996) 693 [arXiv:hep- 
ph/9607310]. L. Covi, E. Roulet, Phys. Lett. B399 (1997) 113 [arXiv:hep-ph/9611425]. A. 
Pilaftsis, Phys. Rev. D56 (1997) 5431 [arXiv:hcp-ph/9707235]. E. Akhmedov, V. Rubakov, 
A. Smirnov, Phys. Rev. Lett. 81 (1998) 1359 [arXiv:hep-ph/9803255]. A. Pilaftsis, T. Un- 
derwood, Nucl. Phys. B692 (2004) 303 [arXiv:hep-ph/0309342]. 

[11] See Y. Burnier, M. Laine, M. Shaposhnikov, JHEP 0602 (2006) 007 [arXiv:hep- 
ph/0511246]. 

[12] The LEP Electroweak Working Group, web page http://lepewwg.web.cern.ch. 



11 



[13] R. Barbieri, P. Creminelli, N. Tetradis, A. Strumia, Nucl. Phys. B575 (2000) 61 [arXiv:hep- 
ph/9911315]. See also E. Nardi, Y. Nir, E. Roulet, J. Racker, JHEP 01 (2006) 164 
[arXiv:hep-ph/0601084]. A. Abada, S. Davidson, F. Josse-Michaux, M. Losada, A. Riotto, 
JHEP 09 (2006) 010 [arXiv:hep-ph/0605281]. 

[14] E. Ma and D. P. Roy, Nucl. Phys. B644 (2002) 290 [arXiv:hep-ph/0206150]. B. Bajc 
and G. Senjanovic, JHEP 0708 (2007) 014 [arXiv:hep-ph/0612029]. B. Bajc, M. Ne- 
mevesek, G. Senjanovic, arXiv:hep-ph/0703080. R. Franceschini, T. Hambye, A. Strumia, 
arXiv:0805.1613. 

[15] M. Miihlleitner, M. Spira, Phys. Rev. D68 (2003) 117701 [arXiv:hep-ph/0305288]. K. Huitu, 
J. Maalampi, A. Pietila and M. Raidal, Nucl. Phys. B487 (1997) 27 [arXiv:hep-ph/9606311]. 
A. Hektor, M. Kadastik, M. Muntel, M. Raidal and L. Rebane, arXiv:0705.1495. E. J. Chun, 
K. Y. Lee and S. C. Park, Phys. Lett. B566 (2003) 142 [arXiv:hep-ph/0304069]. W. Chao, 
Z. G. Si, Z. z. Xing and S. Zhou, arXiv:0804.1265. P. Fileviez Perez, T. Han, G. Y. Huang, 
T. Li and K. Wang, arXiv:0805.3536. 

[16] See e.g. M. S. Carena, A. Megevand, M. Quiros and C. E. M. Wagner, Nucl. Phys. B176 
(2005) 179 [arXiv:hep-ph/0410352] and ref.s therein. 

[17] J.M. Frere, T. Hambye, G. Vertongen, arXiv:0806.0841. 

[18] O tempora, o mores! 

[19] L. E. Ibanez and A. M. Uranga, JHEP 0703 (2007) 052 [arXiv:hep-th/0609213]. 



12 



